Unstable states in QED of strong magnetic fields* 
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We question the use of stable asymptotic scattering states in QED of strong magnetic fields. 
To correctly describe excited Landau states and photons above the pair creation threshold the 
asymptotic fields are chosen as generalized Licht fields. In this way the ofi^-shell behavior of unstable 
particles is automatically taken into account, and the resonant divergences that occur in scattering 
cross sections in the presence of a strong external magnetic field are avoided. While in a limiting 
case the conventional electron propagator with Breit-Wigner form is obtained, in this formalism it 
is also possible to calculate S-matrix elements with external unstable particles. 



I. INTRODUCTION 



The discovery of neutron stars with extremely strong magnetic fields B up to 10^^ Gauss has given the impetus 
to numerous calculations of QED processes in which the magnetic field is taken into account exactly. A remarkable 
feature of these magnetic field strengths is that the cyclotron energy becomes of the order of the electron rest energy, 
r. and, consequently, the quantization of the electron states into discrete Landau levels becomes important. These QED 
^ ■ processes were recalculated using conventional perturbation theory in the Furry picture: the free electron propagator 
0^ ■ and the free wave function of an electron were replaced by the exact Green's function and the exact solution of the Dirac 

■ equation with a homogenous, static magnetic field. Then the same Feynman rules were applied. In this way, most 
' first- and second-order processes have been recalculated in the last fifteen years (for a review cf. Ref. However, 

the results often show an unsatisfactory behavior. Because of strict energy conservation, first-order processes such as 
cyclotron absorption have delta-function-like decay-widths or, like cyclotron emission, become infinite averaging over 

■ "reasonable" distribution functions. This is remedied by accounting for the finite lifetime of the external states; i.e., 
one replaces the delta-function expressing energy conservation by a Lorentz curve in the decay-width. In the case of 

'" p^ 1 second- or higher-order processes singularities arise due to on-shell intermediate states. Here, one accounts for the 
I ' finite lifetime of the intermediate Landau states and replaces in the electron propagator Feynman's ie with ^iF^r^^, 
Qh! where F^v.t is the decay-width of the electron state with Landau quantum number N and polarization r. 
^ ' Despite the use of the electron propagator with a complex mass some processes still lead to divergent cross sections. 

The most prominent example is magnetic Compton scattering with an initial photon which is above the pair creation 
^ • threshold. If the intermediate electron is in the stable Landau ground state, there is no decay-width associated with 
' it, and, consequently, the total cross section is divergent everywhere above the pair creation threshold ||]. Another 
problem arises if more than one particle is unstable. Then it is not obvious how the decay-widths should be defined. 
d Usually, for first-order processes the total decay- width is assumed to be additive, i.e., to be the sum of the decay- widths 
of the individual particles. In contrast, for second-order processes the decay- width of every single virtual particle is 
chosen as its individual on-shell decay-width. This seems to be arbitrary and shows the absence of a comprehensive 
strategy to treat the instability of electrons and photons in magnetic fields. Therefore, it is the purpose of this paper 
to formulate a well-defined perturbation theory for QED of strong magnetic fields where the finite lifetime of excited 
Landau states and photons above the pair creation threshold is automatically incorporated. To this end, we give up 
the concept of stable scattering states and instead introduce generalized Licht fields for the unstable particles. The 
energy of particles described by Licht fields is not fixed by an on-shell condition but is given by some spectral function. 
We do not attempt to calculate these spectral functions from first principles because they are well approximated for 
practical calculations by Lorentz curves ||^. The advantage of using Licht fields is that within this formalism it is 
possible to take into account consistently the instability of intermediate and external particles. 

As an application we show, for the generic example of magnetic Compton scattering in which the initial photon 
is above the pair creation threshold, how the use of Licht fields eliminates resonant divergences of QED with strong 
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magnetic fields. 



II. UNSTABLE STATES AND PROPAGATORS 



To see the underlying reasons for the unsatisfactory behavior of cross sections of QED with strong magnetic fields, 
we remind the reader of two failures of perturbation theory in the Furry picture for B > using stable particle states 
and propagators: 

i) It is a highly distinctive feature of vacuum theory that the Hilbert space structure does not change in going 
from the free to the interacting theory. This remains true for B > only if the electron self-energy and the vacuum 
polarization vanish on-shell. But this is not the case due to the imaginary part of the self-energy indicating the decay 
of Landau states for N > and of the vacuum polarization indicating the decay of a photon in an e^e+-pair above 
the pair creation threshold . 

it) Fields that describe unstable states have a vanishing LSZ-limit 

For us these results are the starting point for identifying the correct states and propagators for decaying electrons 
and photons in strong magnetic fields. There are two ways of introducing unstable particles in quantum field theory. 
Usually, the notation of a complex mass shell, p^p'^ — — iF, is used. This approach is easily applied for propagators 
yielding the typical Breit-Wigner shape for resonances in cross sections but was recently also generalized to 
external lines Another method is the use of generalized Licht fields. Here, one replaces the on-shell condition 
Pf_iP^ — vr? by an off-shell mass spectral density. The abandoning of the on-shell condition is justified by the time- 
energy uncertainty which forbids an unstable particle to have a fixed energy. The use of Licht fields is additionally 
motivated by the following two reasons: First, generalized Licht fields have non-vanishing LSZ-limits. Second, in the 
case of propagators, the Licht field approach is the more general one and contains the complex mass shell method as 
a special case. Therefore, we follow the second approach in this work and, in the spirit of Refs. introduce 
suitable generalized Licht fields for the unstable particles. However, the two different methods result in a different 



treatment of external lines. These differences will be discussed at the end of Section |I] . 
For the electron we define the Licht field by 

/•oc 

and ^ — 5*^7°, where the anticommutation relations are the usual ones for = 0; for > 

[ha{E),h\,{E')\ = [d^{E\d\,{E')\=8,^,.8{E-E') (2) 

and the other anticommutators are zero. Here, V-'a'*''' a-re the energy solutions of the Dirac equation in the presence 
of the external magnetic field B = -Be^, a — {N,T,py,pz} denotes the set of quantum numbers needed in order to 



completely characterize the solutions, and A = ± distinguishes positive and negative energy solutions |11 12]. The 



energy E of the particle is smeared around the on-shell value En — \J rn?- -|- 2NeB + pi due to the integration over 

1/2 

the spectral functions Zy^ ^ (E) . These functions are generalizations of the wave function renormalization constant 
Z of a stable field and labeled by that subset of quantum numbers n — {A, r} which enters in the decay-width of 
the unstable states. In the following, we will suppress the other quantum numbers Py and pz |l^ ]. The dependence 

1 /2 

of Z-^^ on A reflects the different time evolution of positive and negative energy solutions. Therefore, the difference 

1/2 1/2 

between Z_^!^n ^^'^ ^- n should show up only as some kind of boundary condition. In order to get a charge symmetric 

theory, the condition Z^^^"^ = e^'^Z^/t, where is an arbitrary phase, follows. 

A physical interpretation of Eq. (|l]) is that 4" describes n-times different particles, i.e., every Landau state with 
distinct A^ and r would be identified as a different particle. Excited states with A^ > are unstable because of 

1 /2 

the interaction with the photon field. These particles can decay and Z^^ [E] weights the contributions of creating 
an unstable "(A^, T)-particle" with energy E and polarization r. Only the ground state n — {N = 0,t = — 1) 
remains stable. To recover the usual on-shell energy relation and wave function of a stable electron, one has to set 

^A^o^-i ~ ~ ^o)- Then, the ground state has the usual LSZ-limit V'o.-i ^^'^'^a^-i for t oo, where Z2 is 
the normal electron wave function renormalization constant and the limit - as in all of the following - should be 
understood in the weak operator topology. 

In the case of the photon field with its continuous energy spectrum we adopt, in an analogous way, the Licht field 
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/•oo 

A,{x) = / ds' (z]'^,{s')A[+\x, s') + z'J',{s')Ai-Hx, s')) , (3) 



where and A^ ' are the positive and negative energy photon fields, respectively. Similarly to the case of 

1/2 

the electron field, the functions Z-^^ (s') are labeled besides by A by those quantum numbers on which the I-7- 
pair production probability rj.(s) depends: the energy perpendicular to the magnetic field s = tjsin^, and the 
polarization r of the photon . The part of the photon field A^^"* (x) with energy below the pair creation threshold 

1 /2 

uj — 2m/ sinO is stable. Therefore, for s < 2m, the functions Z-^\, (s') = S{s — s') and the field has the usual LSZ-limit 

Ai^\x, s) zl^^Al^^'^^^x, s) for t 00. 

In contrast with the fields 4* and describing unstable particles, the component fields A'^^\x,s') and 

V;W(x,i?) = &„(£;)^(+)(x)e-^* (4) 

^i-Hx,E)=di{E)i;i-\^)e+^^' (5) 

do have a non-vanishing LSZ-limit 

Al^^ (x, s') ^ ix,s') for t^oo (6) 

and 

i>i^Hx,E)^^i^^°-'{x,E) for t-.^. (7) 

Therefore, one is able to compute Green's functions with the Gell-Mann-Low or the LSZ-reduction formula using the 
decomposed out-fields. Expressed in terms of these, the interaction Hamiltonian Hj reads 

d'^Y. J2 dEj dE, ds'zl{l^{E,)zlil^{E,)Zl!l{s') 

where we omit all counterterms. Since here we are only interested in first and second-order processes, we do not 
take care of renormalization. But we want to mention that, since physical parameters should be chosen as directly 
observable quantities, the bare parameter mo should be expressed by the physical mass m only for the Landau ground 
state. In this way, the usual electron wave function renormalization constant Z2 is fixed. However, for decaying states 
one should choose as a physical parameter instead of m some characteristic parameter of an unstable state, e.g. the 
decay- width r„. 

Now, we are ready to derive diagrammatic perturbation theory. As a first step we compute the components of the 
electron propagator (omitting the index "out" from now on) for A'^ > 0, 

iSF{xi,Ei,ni;x2,E2,n2) =< 0|T (^^^^(xi, Si)V^i^^^(a;2, £^2)) |0 >= 

= SiE, - E2)5n,,n. r ^ -j^ ^(^)(xi)^i';(x2) e-(*-*^) , (9) 

27r s - A (£^1 - le) 

and read off the vertex as 

^er4iljE,)Zlil^iE2)Zl{l{s') . (10) 

Using the decomposed propagator, one has not only to integrate over all not fixed momenta, but also over the variables 
of the functions Z. However, we are mainly interested in Green's functions of The total propagator is sandwiched 
between two vertices. One of the two integrations over E breaks down due to the delta-function and one obtains 

iSf{xi,X2) =< 0\T {^ixi)^{x2)) |0 >= 

HdEZ^m r ^ V^^)(x,)?i''(x3)e-(*^-*-) . (11) 

^Jo J^^2n s- X{E -le) 

Here, we set the ill-defined Zq,^i{E) equal to 5{E — Eq) to obtain a compact expression and the vertex becomes the 
usual iej^. The functions Zn{E) are abbreviations for [zj/^l"^. Therefore, they are real and independent from A. One 
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should remember that the propagator obtained, although similar to the spectral representation of the full propagator 
in vacuum theory, is a bare one. 

The derivation of the photon propagator is similar. However, there arises the additional difficulty that the usual 
spin projection operators do not work for off-shell states: a propagator for spin-s-particlcs will generally contain 
particles with lower spin values (s — 1, s — 2, . . . , 0). But since Feynman diagrams with virtual photons do not play a 
prominent role in the astrophysical applications, we omit the derivation here |l5| . 



III. SPECTRAL FUNCTIONS 

In order to make the whole treatment consistent, the spectral functions Zn{E) and Zr{s) have to be - at least in 
principle - computable. A hard way is to use the fact that the Z^'^{E) completely determine through the Eqs. (|l|) 
and (0) the normalization of the components of the field The latter is fixed by the canonical anticommutation 

1/2 

relations. Therefore one can use perturbation theory in the Heisenberg picture to calculate Z-^^{E), and, similarly, 

In practice, these calculations are nearly intractable and one will use a guess. The Ansdtze according to conventional 
wisdom are Lorentzians, 

Zi:{E) - —, ^\ ^ , (12) 



for the electron and 



Z^i^')-—^ T (13) 



for the photon, where we choose r„ and Tr[s) to be the total decay- width of the Landau level {N,t) [|12[ and the 
l-7-pair production probability Tr{s) respectively, calculated in conventional perturbation theory. 

This Ansatz reproduces the electron propagator with the normally used Breit-Wigner prescription ie — > ^ir„ for 
r„ En- In this case, after expanding the poles, the lower limit of integration can be extended from to —00 
producing two poles at s — X{En — ^ir„). According to Ref. the Breit-Wigner approximation yields a result that 
is always consistent within the perturbation theoretical order of calculation. In particular, the author showed that 
for all B and N the deviations from the Breit-Wigner line shape are small. From our derivation follows the usual 
restriction r„ <C _E„ for the validity of the approximation. Using the approximative formulae of Ref. p7[ | for r„, one 
sees that the condition r„ <^ En is indeed always satisfied. 

However, from a more fundamental point of view, it is clear that the extension of the integration to negative E 
violates the spectral condition and is at the root of the violation of unitarity and causality. Although therefore 
the consistency of the approach is lost, it seems to us worthwhile to explore the consequences of this Ansatz. In 
contrast with the usual derivation of the propagator with Breit-Wigner shape, the derivation presented here gives a 
complete scheme that describes unstable particles as external as well as intermediate states. Furthermore, one can 
treat processes with more than one unstable particle without ambiguities. Therefore, in this section we do not attempt 
to calculate the spectral functions Z from first principles but restrict ourselves to the simpler task of investigating 
the consequences of the Breit-Wigner approximation in the Licht field approach. 

Denoting the electron wave functions obtained in this approximation by ipj^liix), we obtain 

i^i^lix) = iV^(^)(x)e-'^^"*-^r"*0(i) (14) 
and ^p'j^l = Here, N is a normalization constant and we chose the signs in 

V n h - En + ^tAL n 

according to the following two requirements: first, we demand that the wave functions do not vanish for t > and, 
second, in the limit r„ ^ the phase of V'^'^i h^is to coincide with the phase of ip^\ The first requirement leads 
automatically to decaying states, for both positive and negative energy solutions. The choice of non-vanishing wave 
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functions for t < results in states whose norm grows in time. Therefore, the choice between non-vanishing wave 
functions for t > or i < corresponds to the choice of the direction of the time arrow and has to be made by hand. 

1/2 

As anticipated, the different form of „ for A = ± is necessary to obtain the correct boundary condition for decaying 
states. Similarly, we obtain for the photon wave functions 

andA(-;U4+;)*. 

Since the norm of the states varies with time, the correct normalization is not obvious. A reasonable prescription 
is the requirement that cross sections calculated with decaying states coincide in the limit F ^ with the same cross 
sections calculated in the usual formalism. 

Formally, S'-matrix elements with external unstable particles will be calculated in the usual way. But because of 
the 6 functions, the time integration over vertices with external unstable particle goes effectively only from to oo. 
Thereby, no divergent time integrals will be caused by the real part of the exponentials of decaying states. 

Finally, we want to compare our approach with that of a recent paper |^ . The authors propose the use of propagators 
with complexified energies and - this is the main difference - the use of external states 

C(^)=V^i'Hx)e+'^^"*+*^^"* (18) 

for the electrons and similar ones for photons. In contrast to our wave functions, the norm of the negative energy 
solutions grows in time. Furthermore, these wave functions are valid for t < as well as for t > 0. Therefore, 
the real part of the exponentials will lead to divergent time integrals. Consequently, no scattering amplitudes in 
the normal sense (i.e. for transitions between tij ±oo) with external unstable states can be calculated in this 
approach. Instead, one calculates matrix elements of the time evolution operator U{tf,ti) which are dependent on 
the time lapse tf — ti between preparation of the initial state and measurement of the final state. The authors of 
Ref. claim that this approach, which is the direct transcription of the Wigner-Weisskopf method of non-relativistic 
quantum mechanics, eliminates all divergences of QED with strong magnetic fields. However, the main object of field 
theory, the S'-matrix, is not generally computable in this formalism. Moreover, the wave function '0 is not the Dirac 
conjugate spinor of ip since the real part of the exponent changes sign. Consequently, charge symmetry is lost and 
the S'-matrix elements lack crossing symmetry. Finally, a more practical objection seems to be important. Since the 
main application of QED of strong magnetic field is astrophysics, the usefulness of cross sections which are dependent 
on the time lapse between "preparation" and "measurement" of the states is restricted. 



IV. APPLICATIONS 



In this section, we want to illustrate some basic consequences of this formalism. 
First, we consider a generic first-order process. The S'-matrix element is given by 

S/, = J dE^dE2ds' zl'\E^)Zl'\E2)Zl'\s')Sfr , (19) 

where SJf" is the conventional S-matrix element but with off-shell energies. Choosing the spectral functions to be 
Lorentz curves, there is no difference between the approach presented in this work and conventional perturbation 
theory where the delta function expressing energy conservation is replaced ad-hoc by a Lorentzian. Since the wave 
functions depend on F only through an exponential factor, the simple assumption that the total decay-width is 
additive, is valid, i.e. (as used e.g. in Ref. [p^), 

rtot = r„, +r„, + r,(s). (20) 

Hence, quantum correlations between different decaying states exist only for non-Lorentzian spectral functions. 

Second, let us consider magnetic Compton scattering as a typical second-order process. The electron propagator 
coincides with the conventional one with complexified energy for Zn{E) = Z^{E). Therefore, in the simplest case 
of A'i = ^ Nf = Compton scattering where the energy of the photons is below 2m/ sin 0, we obtain the old, 
well-known result Otherwise, the instability of the initial and final particles is also incorporated in the S-matrix 

element. Assuming the functions Z^/^ are well-behaved, the integration over the off-shell energies will remove the 
remaining singularity of the electron propagator when the virtual electron is in the stable Landau ground state. In 
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particular, the S-matrix element of magnetic Compton scattering is now finite in the case where the initial photon is 
above the pair creation threshold. 

Now we want to make our argument quantitative. The 5-matrix element of A^i = ^ Nf = Compton scattering 
is given by 

Sf, = {ief j d^xd^x'4,f\x)-t,,iSF{x,x')-(,^\'^\x') j ds[ds)z]')^{s[)z]!^^^{s'f) 

where i and / refer to initial and final states, while the the quantum numbers of the virtual electron will be marked by 
the subscript a. Using stable fields, the divergence for Si > 2m occurs in the exchange diagram 5^^^ when the virtual 
electron propagates as a positron in the Landau ground state, i.e. has the quantum numbers Na = and Aq = — 1 
Since the space integrals remain unchanged and are finite, we only have to consider the time integrals of 5^^-' for 
TVq = and Xa = 0, 

5(2) oc / dtdt'e'^^' (e^E^it-t-) f il ^^^t^t')] ^-rEa' 



27r 

ds:d44(,2^(s:)zi/^^.(4)e— '*'e-/* = 5f ) . (22) 

Here, E denotes the energy of the electrons and lo' = s' sinO the (off-shell) energy of the photons. Performing the two 
time integrations and the integration over 2, which comes from the ^-function of the electron propagator, results in 

5p) ^ / ds^ds', ^^i:!^^^ 2nSiE, +.',-E.- . (23) 
The result of the integration over s'^ depends on the energy perpendicular to the magnetic field of the final photon. 



.(2) _ 2^4(n(gb) 



Sy> = ' sin^ 



-Ei - Ea + ojf + ie 



Sf < 2m 



^■'^ -Ei - Ea + co'j: + is 



s'j > 2m ' 



where we introduce s'^ i = (Ef + ujf — Ei) sin 9i, 2 — (^/ + '^z ~" -^0 sin^i and loj denotes the on-shell energy of 
the final photon. In the second term of S^^\ the integration over s'f gives a finite result as long as a principal value 
integral of the integrand can be defined. Therefore, the only remaining dangerous part of S^^'^ is the first term. But 



for s/ < 2m the denominator can never become zero because Ei+ Ea — m + y m^ + cos^ 0/ is always greater than 

ujf. Here, we assumed pi_z = without loss of generality and used momentum conservation parallel to the magnetic 
field. 

In the more general case of Compton scattering when Ni and Nf are not restricted to be zero, this result remains 
valid. In this case, there are additionally integrations over the off-shell energies of the unstable electrons. As above, 
the necessary condition for a finite ^-matrix element is that a principal value integral of the integrand can be defined. 

This example illustrates well the connection between the instability of external particles and singularities of inter- 
mediate states: as soon as the energy of the initial photon is above the pair creation threshold, the virtual electron 
can become real, producing a divergent cross section. At the same time, however, the photons also become unstable. 
Taking this instability into account properly, one obtains well-behaved cross sections. 

Finally, we want to comment on the behavior of magnetic Compton scattering if all external particles are stable. 
Then there is no decay-width which could cure the resonance if the virtual electron is in the Landau ground state 
Na = 0. Formally, the resonance energies are given for all Na = 0, 1, ... by pi 



^rcs 



(m^ -1- 2NaeB sin^ 6,) ^'"^ - ml / sin^ , (25) 



i.e. i^res^° — 0- the case of iV^ ^ ^ iV/ — Compton scattering, the limit ^ results in divergent cross 
sections In Ref. |2^, these divergences were interpreted not as resonances but as infrared divergences. By 

contrast, in the case of A'i = ^ N f = scattering, the S'-matrix element diverges like w"^ while the cross section 
goes to a finite, constant value in the limit cj^ — > 
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V. SUMMARY 



We have presented a consistent method to describe the instabihty of excited Landau states and photons above the 
pair creation threshold in QED of strong magnetic fields. This approach consists in using Licht fields for unstable 
states and introduces additionally integrations over the off-shell energies of the unstable particles. We have shown for 
the generic example of Compton scattering where the energy of the photon is above the pair creation threshold how 
in this way the resonant divergences of S'-matrix elements of QED of strong magnetic field are avoided. 

In the Breit-Wigner approximation the Licht states are exponentially decaying or growing in time. Since the 
divergent part of the wave functions is cut off by 9 functions, no divergent time integrals will be caused by the real 
part of the exponentials. Therefore, in this formalism it is possible to calculate S'-matrix elements with external 
unstable particles. 
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